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ABSTRACT 

We present the results of a deep photometric study of the outer halo of NGC 
1275, the highly active cD galaxy at the center of the Perseus cluster. We find 
a modest excess of faint {R > 22.5) starhke objects in its halo, indicating a 
population of old-halo globular clusters. However, the total estimated cluster 
population corresponds to a specific frequency of = 4.4 ± 1.2, no larger 
than that of normal giant ellipticals and three times lower than that of other 
central cD galaxies such as M87. We discuss several ideas for the origin of this 
galaxy. Our results reinforce the view that high Sn (ie: highly efficient globular 
cluster formation) is not associated with cooling flows, or with recent starburst 
or merger phenomena. 
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1. Introduction 

Galaxies are known to exhibit wide differences in tlie relative numbers of their globular 
clusters (GCs). A useful quantitative measure of this effect is the specific frequency (Harris 
& van den Bergh 1981; Harris 1991), 

SN = Nf 10°-^(^^+^v) (1) 

defined as the ratio of total cluster population Nt to galaxy luminosity My. Sn represents 
the overall efficiency of cluster formation (or accretion) averaged over the galaxy's history. 
Spiral galaxies tend to have the smallest populations (^jv ~ 1 — 3), while normal ellipticals 
rank somewhat higher {Sn ~ 4 — 6; see Harris 1991). In this respect, an outstandingly 
interesting group consists of the cD galaxies found in the centers of rich clusters. Although 
some of these have Sn ratios comparable to ordinary ellipticals, others have ratios as high 
as Sn ~ 15. The suggested origins for these very large GC populations have included 
formation out of cooling fiows (Fabian, Nulsen, & Canizares 1984), mergers of spirals 
(Schweizer 1987; Ashman & Zepf 1992; Zepf & Ashman 1993), accretion by tidal stripping 
and cannibalism of neighboring galaxies (e.g. Muzzio 1988), or a high GC formation 
efficiency during the protogalactic era (Harris 1991). Each of these options has distinct 
difficulties in explaining the data: a thorough summary of the issues for these cD galaxies 
is given by Harris, Pritchet, & McClure (1995 = HPM95). 

NGC 1275 (Perseus A) potentially provides an interesting touchstone for many of 
the arguments surrounding the issue of GC formation, since it possesses a combination of 
properties that are rarely found even among central giant ellipticals. It is an extremely 
luminous cD-type galaxy at the center of the rich Perseus cluster (Abell 426), it has one 
of the largest known cooling flows (~ 200 M0/yr; Fabian et al. 1981), a highly active 
Seyfert-type nucleus, and has undergone very obvious GC formation in its central regions 
in the past ~ lO'^ y (Holtzman et al. 1992; Richer et al. 1993). Zepf et al. (1995) note 
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the unusual distended shape of this Seyfert galaxy and suggest recent merger events as 
the cause. Many other features add to the interest of this galaxy. A system of luminous 
filaments (see Fig. 0) stretches more than 20 kpc from the galaxy's center. These filaments 
have a systemic radial velocity of ~ 5300 km s~^, the same as the main body of NGC 1275 
(the 'low velocity' or LV system). They may have been born out of the cooling flow and 
ionized either by nonthermal emission from the nucleus (Kent & Sargent 1979) or by shocks 
within the LV material itself (see Heckman et al. 1989 for a review). The nucleus is a region 
of high Ha emission and is host to both molecular (Lazareff et al. 1989) and neutral (Jaffe 
1990) hydrogen. 

Equally important is the existence of a 'high velocity' (HV) system of giant HIT 
clouds with Vo ~ 8300 km s^^ which occupies the N and NW regions of NGC 1275. 
Various theories exist regarding the relationship of the LV and HV components. The 
initial supposition (Baade & Minkowski 1954) was that the HV system represents the 
gaseous remains of another galaxy, perhaps a large spiral, merging with the underlying 
giant elliptical. Alternately, De Young, Roberts, & Saslaw (1973) proposed that the two 
components constitute a chance superposition of two galaxies that do not interact. Still 
another proposal (Burbidge & Burbidge 1965) is that the HV system was ejected from the 
core of NGC 1275. 

Later evidence has served to place the HV system at least partially in the foreground: 
HI associated with the HV system has been found in absorption in the spectrum of the 
nuclear radio source 3C 84 (De Young, Roberts, & Saslaw 1973), and IlJE spectra show 
HV-related Lya in front of the nucleus (Rubin et al. 1977; Briggs, Snijders, & Boksenberg 
1982). Similarly, absorption of NGC 1275 starlight has been associated with the areas of 
strongest HV emission (N0rgaard-Nielsen et al. 1993). This evidence, combined with the 
lack of a symmetrical counterjet has made the ejection hypothesis unlikely. New evidence 
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for the interacting-galaxy hypothesis is provided by Boroson's (1990) discovery of stellar 
absorption in the Ca II A8948 line of the HV system, together with a new velocity curve 
that suggests a velocity gradient across the HV gas. 

The recent discovery of several compact blue objects in the galaxy core (Holtzman et 
al. 1992) has attracted much attention. The initial identification of many of these objects 
as young (< 1 Gyr) globular clusters has been substantiated by later photometry and 
spectroscopy (Richer et al. 1993; Zepf et al. 1995). Condensation out of the cooling flow, 
past mergers, and encounters with the HV system have all been suggested as possible causes 
for their origin. Holtzman et al. (1992) and Zepf et al. (1995) contend that the young GCs 
are evidence of a merger with another galaxy, which would have brought in a fresh supply 
of relatively cool gas. It is not obvious whether or not the HV system represents part of 
the merger remnant, although Richer et al. (1993) point out that the high relative velocity 
of the two systems (~ 3000 km/s) gives an interaction timescale of only ~ 10^ yr, which 
does not agree with the proposed cluster ages of 0.5 Gyr (Zepf et al. 1995). In addition, 
examination of the velocity curve for the HV system suggests that it is only just beginning 
to come in contact with the main body of NGC 1275 (Boroson 1990). Lester et al. (1995) 
identify thermal emission from a dust component and claim that its presence argues in 
favor of a recent merger. Thus, the nuclear activity may be the visible result of a merger 
that took place in the last ~ 1 Gyr and is by now nearly complete, with the HV system 
being a still more recent arrival. 

//globular clusters form efficiently out of spiral mergers (e.g. Zepf & Ashman 1993) or 
coohng flows (Fabian et al. 1984), and if the very rich cluster systems in cD galaxies built 
up this way, then we might anticipate that the 'old', outer halo of NGC 1275 would contain 
a populous GC system. An additional reason for expecting NGC 1275 a priori to have a 
large cluster population is the strong correlation between high and environment. Of the 
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half-dozen known high-S'jv systems (HPM95; Cote 1995), all are very near the dynamical 
centers of rich galaxy clusters. Specifically, these are the cD galaxies in Virgo, Fornax, 
Hydra I, Coma, A2052, and A2107. NGC 1275 is an extremely luminous galaxy quite near 
the dynamical center of the Perseus cluster: at a redshift Vo — 5316 ± 78 km/s, its velocity 
is little different from the mean cluster redshift of 5451 ± 168 km/s (Bird 1994). To compare 
this peculiar velocity to that of similar cD galaxies. Bird uses Gebhardt and Beers' (1991) 
Z score parameter 

ZoB = (2) 

where Cbi and Sbi are the biweight estimators of mean cluster velocity and velocity 
dispersion (these are more robust than the usual Gaussian centroid and dispersion, which 
tend to fail when applied to non-Gaussian distributions). Bird lists the Z-score for NGC 
1275 as Z = —0.108 (corrected for substructure), i.e. roughly a tenth of a standard 
deviation from the Perseus centroid. When we add the fact that NGC 1275 is at the center 
of an extensive X-ray halo, the evidence is strong that it is sitting very nearly at the center 
of Perseus. By association with the other cD galaxies mentioned above, its old-halo cluster 
population should then be very rich. If, on the other hand, its specific frequency is low, 
then several other models for GC formation in such galaxies might be viable. 

There is little information available on the old-halo clusters in NGC 1275: to date, 
the sole study addressing the old GC population is that of N0rgaard- Nielsen et al. (1993). 
Prom images of the inner 1' of the galaxy, they suggest by extrapolation that NGC 1275 
may have a total GC population 2/3 as great as that of M87 (Nt ~ 10^), in which case it 
would have a normal Sn- However, the internal uncertainties in their data are too great to 
permit definitive conclusions. 

In this paper, we present the results of a deeper and more extensive survey of the outer 
halo of NGC 1275, and discuss its imphcations for GC formation scenarios. 
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2. Observations and Data Reduction 

The images for this study were obtained 1993 September 21 at the prime focus of the 
Canada- France-Hawaii Telescope (CFHT), with the High Resolution Camera (HRCam; 
McClure et al. 1989) and the Loral3 CCD binned 2 x 2 to give an image scale 0.216 
arcsec pixel"^. The usable photometric field of the camera was 750 x 750 binned pixels, or 
2.7 arcmin on a side. A 20 x 300 sec exposure series in R nearly centered on NGC 1275 
was obtained as the primary material for our study (attempts to obtain comparably deep 
images in other filters, as well as a background control field, were unfortunately terminated 
by weather and instrument problems). The R band was chosen because it is near the 
peak sensitivity of the CCD, and is well matched to the moderately red colors of our 
target objects (the old globular clusters). The series of images was re-registered and 
median-combined to yield the composite picture shown in Fig. Q the mean FWHM (seeing) 
on this combined frame is 3.4 px (0.75 arcsec). 

Photometric calibration was carried out by transferring short exposures of Landolt 
(1992) standard stars onto our individual NGC 1275 exposures. The resulting solution for 
the standard stars was 

R = r- O.lOX - 0.013(y - R) ^ const (3) 

where r is the instrumental magnitude through a large (15 px) aperture and X is the 
airmass. This relation reproduced the standard-star values to within ±0.01 mag. Since 
the color coefficient is very small and we have only one bandpass for our program field, we 
neglected the V — R term. Final PSF-fitting photometry of all the objects on the composite 
frame was carried out only after removal of the overall contours of the galaxy light. A model 
for the isophotal contours was created using the STSDAS routines isophote. ellipse and 
isophote.bmodel. This model was then subtracted from the picture to eliminate the major 
part of the light gradient from the body of NGC 1275. A sequence similar to that described 
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by Fischer et al. (1990) was then implemented, whereby passes of DAOPHOT/ALLSTAR 
(Stetson 1992; Stetson, Davis, & Crabtree 1990) were alternated with smoothing and 
removal of the remaining background light with the fast ring- median filter of Seeker (1995). 
Two iterations of this sequence proved to be sufficient to locate virtually all of the starlike 
objects in the field. 



3. Analysis 

To evaluate the photometric completeness as a function of magnitude and position on 
the frame and to help set parameters for rejection of nonstellar images, we carried out a 
series of simulations. First, all the detected objects and the galaxy light were subtracted 
to produce a frame bare of any features except background sky noise. With this blank 
frame we then recreated a dozen simulated frames, adding a total of 3600 stars in groups of 
300 over a wide range of magnitudes. These were spatially distributed to mimic the light 
profile of the galaxy (see below), following a simple power-law distribution centered on the 
main galaxy. The results of these experiments for the completeness fraction / (the fraction 
of objects recovered at any magnitude) are summarized in Figure ^ and Table 1. Fig. |^ 
displays the experimental results, fitted by Pritchet's interpolation function (see Fleming et 
al. 1995), 

^l + a^R-Ri,my\ 

The limiting magnitude (i.e. the level at which / drops to 50%) is i? ~ 24.43, a value 
nearly independent of radial location for our main region of interest (r > 100 px, outside 
the central region where the galaxy light dominates; see below). All fainter objects were 
discarded from our lists. Note that the completeness function drops quite steeply, so that 
for R < 24 the photometry is virtually 100% complete. 
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For any galaxy at the distance of NGC 1275, globular clusters appear as faint starlike 
objects sprinkled around its halo. Yet due to the presence of the optical filaments, the 
NGC 1275 field is photometrically much more complex than that of a more normal giant 
elliptical. These features caused the DAOPHOT algorithms to fit PSFs to the clumpier 
regions of the filaments, thus generating numerous spurious detections. When these are 
added to other types of contaminating objects on the frame (faint background galaxies, 
Perseus dwarfs, etc.), the majority of the detected objects in the raw lists turn out to be 
nonstellar and need to be eliminated as completely as possible. To accomplish this, we used 
the radial moments ri and r_2, which measure image wing spread and central concentration 
respectively, along with an isophotal magnitude Ri (Harris et al. 1991). These plots 
are shown in Figure ^ where the starlike objects lie along narrow sequences. The three 
parameters are so well correlated in this field that they are largely redundant, and the 
criterion we adopted was simply to reject objects with ri > 2.5, since this moment gave the 
clearest empirical discrimination of nonstarlike objects. This step eliminated most of the 
extended objects and multiple detections from the list. Finally, the remaining 349 objects 
(candidate GCs) were then examined one by one against the original image. All objects 
within 5 pixels of the frame edge, or which were associated with bad pixels, filaments, or 
were otherwise anomalous to careful inspection, were also rejected. This culling procedure 
left 294 final candidates. 

Since the old-halo globular clusters have such a similar luminosity function (GCLF) 
among all types of large galaxies (especially among the giant ellipticals; e.g., Harris 
1991; Harris et al. 1991; Ajhar et al. 1994; Whitmore et al. 1995), we can make an 
excellent preliminary guess at the magnitude range over which we expect to detect a 
cluster population around NGC 1275. For the well known Virgo ellipticals the GCLF 
peaks at i? ~ 23.2 ± 0.2 (or A/r = —7.8 for (m — M)y = 31.0) and with a dispersion near 
(T = 1.4 ± 0.1 (Harris et al. 1991; Whitmore et al. 1995; Harris, Harris, & McLaughlin 
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1996 Seeker & Harris, 1993). At a redshift Vq ~ 5400 km/s, NGC 1275 is four times farther 
away; after addition of foreground absorption {Eb-v = 0.15; Burstein & Heiles 1984), we 
would then expect the GCLF turnover to he at i? ~ 26.5, more than two magnitudes fainter 
than our photometric hmit. The GCLF should then start appearing only for R > 22.5 and 
should rise rapidly faintward of that level. Our directly observed result (see below) is fully 
in line with this expected range. 

Figure ^ shows the positions of the final 294 GC candidates over the target magnitude 
range 22 < R < 24.43. Several gaps in the spatial distribution correspond with the 
positions of the bright foreground stars, other large Perseus galaxies, and the HV region 
just northwest of the galaxy core. The latter has been associated with regions of high 
absorption (N0rgaard-Nielsen et al. 1993). 

The elliptical annuli in Fig. | were used to calculate the radial distribution, excluding 
the inner two annuli, for which the number statistics were poor. This distribution is 
listed in Table 2 and displayed in Figure |^, where the mean number density a of the 
detected candidates (corrected for areal and magnitude incompleteness) is plotted against 
radius. A modest falloff in a toward larger radii is evident, suggesting that we have 
successfully detected the presence of an old-halo globular cluster system. However, the 
fact that this system was not obvious from casual inspection of the images (or indeed from 
inspection of the 'cleaned' distribution in Fig. ^), already suggests that NGC 1275 is not a 
high-specific- frequency galaxy like M87. This suspicion will be borne out in the calculation 
of to be made below. 

Next we derive the luminosity function. The measured GCLF is shown in Figure 
^. Here, the outermost three annuli are again used to set the background LF, shown 
by the dashed line. As expected, the residual population (j){R) (number of clusters per 
0.2-magnitude bin) is negligible for R < 22.5 and rises steadily and steeply to fainter 
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magnitudes as is typical for a normal old-halo GCLF. Although the faintest bin shows a 
rather sudden decrease in 0, this does not represent the turnover point, which would be 
~ 2 mag fainter. In this last bin, the completeness fraction is near / ~ 0.5 and is clearly 
the least certain point, although it is still puzzlingly low. Adopting the turnover level and 
standard deviation of the GCLF with uncertainties as described above, we find that the 
observed population brighter than i?(lim) = 24.4 represents only the top (3.75 ± 1.00)% of 
the total GCLF population. 

Since we had no background comparison field, a preliminary estimation of the 
background level was obtained simply by taking the total number of objects in the 
outermost three annuli and dividing by the total area of these, giving (Jb = 39± 11 arcmin"^. 
An alternate approach (e.g. Harris 1986) is to fit a simple power law 

(7 = + ab (5) 

to the data points and to solve for as many of the parameters as possible directly from the 
fit. Since the radial region under scrutiny is not extensive (0.5'< r < 2.0'), the data do not 
constrain the exponent a in any useful way. Instead, we set a by using the relation (Harris 
1993) between a and the integrated absolute magnitude of the galaxy My: 

a = (-0.29 ± 0.03)M^ - 8.00. (6) 

[NB: The zero point of this relation appears incorrectly in Harris (1993). Here we quote the 
corrected version.] 

The integrated absolute magnitude of NGC 1275, with Vt — 11.59 (de Vaucouleurs, de 
Vaucouleurs, & Corwin 1976), is My — —23.3 for a distance modulus (m — M)y = 34.9 
{Hq = 75 km s^^ Mpc~^ and Ay = 0.6). However, the spectrum of NGC 1275 is strongly 
affected by blue light from starburst activity in the core, rendering the galaxy brighter than 
a similar gE in which no star formation was occurring. Romanishin (1987) concludes that 
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the ratio of excess blue light to that of the old halo is 15%. Correcting for this, we will 
adopt My{corr) = —23.14, which produces a ~ —1.3. 

Since the scatter around the mean relation is ±0.2 (Harris 1993), we then fit equation 
(P) using a = (—1.1, —1.3, —1.5). Assuming a, we then solve for and /9 by straightforward 
least squares. In what follows, we use unweighted solutions. Weighted solutions would 
slightly lower ab and leave the effect of increasing the residual cluster populations by about 
20%. 

Fig. 1^ shows all three solutions and the radial profile. The value (Tb ~ 33 arcmin"^ (for 
a = —1.3) may still be an overestimate of the true background level (2' corresponds to a 
linear distance of ~ 40 kpc, and the halo GCS can extend detectably beyond that radius 
for particularly large systems, such as M87). However, lacking a direct determination from 
an adjacent control field, we cannot isolate ab any more narrowly. Statistical estimates of 
CTfe gathered by modelling starcounts and background galaxies are not of any additional help 
because of the unusual complexity of the Perseus field in the neighbourhood of NGC 1275. 

Next are the corrections to the total population due to area. To obtain the total 
cluster population over all radii we integrate relation ^ from 0.5' out to a generous limit 
5' (corresponding to 100 kpc), multiply by the ellipticity area factor q = 0.8, and subtract 
the background value. This calculation yields the values found under Nint of Table 3. A 
rough estimate of the population in the innermost region r < 0.5' may be found by taking 
a constant (background-subtracted) value of cr ~ 40 arcmin^^, which gives an additional 
~ 30 clusters brighter than our limiting magnitude (we emphasize again that this inner 
population is the old-halo GCS, and does not include any of the blue objects discovered in 
the HST survey). Adding these 30 clusters and dividing by (0.0375 ± 0.010) to account for 
the portion of the GLCF which is fainter than our limiting magnitude, we then derive the 
total old-halo cluster populations given in Table 3. 
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For the specific frequency we tlien obtain Sn = 4.3 ± 1.4 - witliin tlie range of 
the Virgo and Fornax elhpticals (Harris 1991), but more than 3 times smaller than the 
'high-^Af' cD systems of which M87 is the prototype. Perhaps a more dramatic illustration 
of this conclusion is shown in Figure in which we show the observed radial profile of the 
M87 GCS as it would appear at the distance of NGC 1275, normalized to the same total 
luminosity. If, as noted above, we have significantly overestimated the background level 
ab, the results for the total population would change, but not dramatically so. Even our 
shallowest adopted profile (a = —1.1) excludes NGC 1275 from the high-S'iv ellipticals. 
Unless deeper and wider-field data reveal major surprises, it seems clear that as giant 
ellipticals go, NGC 1275 has a normal or even subnormal specific frequency. 



4. Discussion 

All previous observations have shown that the conditions necessary to form a high-S'Ar 
galaxy are dominant size and central location within a rich cluster (HPM95). On these 
grounds, we might reasonably expect NGC 1275 to be such a candidate. In addition, 
McLaughlin, Harris, & Hanes (1994) conclude that there is a moderate correlation between 
Bautz-Morgan class and Sat. The BM class of Perseus (II-III) would then also suggest that 
its central cD should be well supplied with GCs. Furthermore, NGC 1275 is undergoing 
highly active phases of gas accretion and star formation. Despite all these conditions, 
its specific frequency is unexceptional. This discrepancy casts some doubt on the view 
(Ashman & Zepf 1992, Zepf & Ashman 1993, Zepf et al. 1995) that merger or starburst 
systems are particularly favorable places for GC formation. 

However, our result strongly reinforces the conclusion of HPM95 that cooling flows do 
not produce large numbers of globular clusters in cD galaxies. In this one respect, NGC 
1275 is similar to certain other cD systems that they discuss. By inference, it would then 
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seem more likely that its very extensive nuclear activity and young cluster formation are 
due to accretion or infall of gas from another source. 

To explain our result for NGC 1275, we are left with a variety of options: 

(1) NGC 1275 was a normal galaxy that has undergone several mergers. In this 
interpretation, we would assume that NGC 1275 started as an unexceptional galaxy (either 
spiral or elliptical) which by chance was close to the dynamical center of A 426. Since 
then, it has merged with numerous other galaxies, producing an overall cD-like and highly 
luminous product in which the traces of the original galaxy have been all but erased. The 
merger or accretion currently taking place is then only the latest in a series. 

Spiral galaxies have low specific frequencies (S'jv ~ 1 — 3; Harris 1991) and merger 
products of spirals are likely to be similarly endowed. It is unquestionably true that 
globular clusters do form in such mergers (Schweizer 1987; Ashman & Zepf 1992; Whitmore 
& Schweizer 1995). But this, by itself, will not necessarily increase the specific frequency of 
the product galaxy. The reason is essentially that ^'at is a population ratio of clusters to field 
stars; therefore, any increase in Sat requires the preferential iorm.&i\OTi. of GCs over field stars 
(which also form during the merger in large numbers of unbound clumps and associations). 
Thus, the final Sn depends critically on how efficiently the merger converts the incoming 
gas into clusters, as well as on the specific frequencies of the incoming pre-merger galaxies. 
The cluster formation efficiency is likely to depend on the detailed geometry of the merger, 
the relative velocities and orientations of the incoming galaxies, the amounts of gas within 
each, and other factors. 

Averaged over many encounters, it may be reasonable to expect a resulting S'jv in the 
range we observe for NGC 1275 only if the formation of GCs was remarkably efficient and 
extremely large amounts of gas were brought in by the mergers (Harris 1995). As yet, no 
models exist to calculate these efficiencies. Although it may not be ruled out, we regard 
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this option as an unlikely solution. 

(2) NGC 1275 was a high-SN giant elliptical galaxy that has undergone many accretions. 
This idea is essentially a more conservative version of the previous one. Here, the initial 
gE might have had a normal specific frequency to go along with its centrally dominant 
position. But since any accreted galaxies (especially spirals) would have had lower specific 
frequencies, the eventual product cD would have been left with a more modest cluster 
population and an S'jv value that has been diluted down through the mergers by a factor of 
2 or more. This idea simply would permit a much lower GC formation efficiency. Again, the 
merger we are seeing today would be only the latest accretion event in an ongoing series. 

While both of these two scenarios arc qualitatively consistent with the observations, a 
few simple calculations cast some doubt on them. If the product galaxy resulted almost 
totally from spiral mergers, then the present-day galaxy with its ~ 8000 old-halo clusters 
would have required the merger of ~45 large spirals the size of the Milky Way (each 
with 100 — 200 GCs) if no new globulars were produced. This is an uncomfortably large 
requirement. However, it is more likely that some GC formation occurred during the 
mergers. If we suppose that roughly half the GCs were formed out of incoming gas in this 
manner, then we might need ~ 25 such spirals, each bringing in ~ 2 x IO^Mq of gas for a 
reasonable GC formation efficiency (see Harris 1995). These numbers translate to a rate 
averaging at least 3 mergers/ Gyr. While this is not impossible, one must assume a very rich 
initial population of big spirals and a central galaxy adept at capturing them. In addition, 
recent evidence (Andreon 1994) shows that previous cluster surveys have underestimated 
the number of spirals by a factor two to three. If the Perseus cluster presently contains 
many spiral galaxies, and the above scenario is correct, the initial population must have 
been overwhelmingly spiral-rich. 

Conversely, if the original central object was a large elliptical containing half the 
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present GC population, then the necessary spiral accretion rate could be reduced to a less 
demanding 1/Gyr. Yet to possess an original specific frequency ~ 15, proto-NGC 1275 
would also be required to be > 4 times less luminous than it is today. If GC formation 
during mergers was significant, then the original galaxy must have been even smaller, 
and alternative (2) then becomes almost indistinguishable from (1). In short, changing a 
high-S'jv galaxy to a Iow-Sn one puts rather severe demands on any accretion or merger 
process. 

(3) NGC 1275 began as a normal-S]^ cD or giant elliptical. This scenario is the most 
conservative of all, since it reduces the need to invoke mergers at any high rate. Rather than 
having accretions dilute or enhance the S'at, we suggest that NGC 1275 began as a gE with 
a specific frequency similar to that observed today. The mergers required to build up the cD 
envelope need not have altered the dramatically, as might be the case if other gas-poor 
gEs of similar Sn were swallowed by the proto-NGC 1275. The present-day activity (nuclear 
star formation, the HV system, etc.) may then be interpreted as fairly rare occurrences 
which will have little effect in the long run. The young blue clusters of Holtzman et al. 
(1992) may dissolve away in the next ~ 1 Gyr from tidal shocking, dynamical friction, and 
self-disruption before they can add to the old-halo cluster population. 

Evidence which could verify or refute the above possibilities may be found in the color 
distribution of the GCS. Unimodality would indicate single- burst formation as in the case of 
ELS-model formation (Eggen, Lynden-Bell & Sandage 1962), while bimodality could point 
to multiple-epoch formation. Furthermore, if a large fraction of NGC 1275's GCs came 
from spirals originally, then they should be metal-poor and thus show up as a significant 
population of bluer objects than normal for a giant elliptical. 

None of the alternatives discussed above is entirely satisfactory, since the small old-halo 
cluster population in NGC 1275 remains unusual in the context of other cD-type galaxies. 
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Although central cDs with normal or subnormal values (< 5) are known elsewhere, none 
have the level of central activity that NGC 1275 exhibits (see HPM95). For further clues, it 
is natural to turn to the combination of features associated with the central regions of this 
object. It is agreed to be a Seyfert galaxy (Vermeulen et al. 1994, Nesterov et al. 1995) 
with a blazar-like spectrum (Lester et al. 1995). Its highly active nuclear region contains 
the radio source 3C 84 which, after a burst of activity observed in 1980, was the brightest 
extragalactic source at mm wavelengths. Vermeulen et al. (1994) have been observing 
the nucleus in the 22 GHz range since 1981 and have recently discovered a counterjet to 
complement the previously known relativistic nuclear jet. Vermeulen et al. believe that the 
counterjet eluded discovery for some time due to a toroid of gas about the central engine 
(assumed to be a black hole) that obscured the counterjet through free-free absorption. 
Presently, the core luminosity is ~ 10^^ ergs s~^, a value in the same range as more distant 
radio galaxies. 

Recent theoretical modelling (Harris & Pudritz 1994) makes the case that globular 
clusters form out of cool gas collected within supergiant molecular clouds (SGMCs). An 
obvious speculation is then that if the nucleus of NGC 1275 was even more active at an 
earlier epoch, the high ionizing radiation and radiation pressure could have disrupted the 
normal GC formation process. That is, by the time the nuclear activity died down to the 
point that star formation could proceed more vigorously, the host SGMCs might already 
have broken down into smaller units within which the formation of massive star clusters 
was impossible (see Harris & Pudritz). However, the problem with invoking this sort of 
process to inhibit GC formation is that most (and probably all) of the central giant cD's 
should have central massive black holes (e.g. Fabian & Rees 1995; Harms et al. 1994; Ford 
et al. 1994). That is, most or all of them should have experienced highly active early phases 



separate phenomenon from that discussed in the Introduction. 
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of disruptive nuclear activity. Yet some (M87, for example) successfully produced huge 
numbers of GCs, while others (such as NGC 1275) could not. In addition, it is plainly true 
that massive star clusters have formed quite recently in the inner few kpc of NGC 1275 
despite its strong nuclear activity; thus it is not clear just how these seemingly disruptive 
effects actually influence GC formation. 

Another piece of the puzzle may lie in the details of the Perseus environment and the 
X-ray halo gas, which provides a signature of the evolutionary state of the cluster. Actively 
studied at X-ray and optical wavelengths, A 426 was at first assumed to be relaxed (Jones 
& Forman 1984). However, more extensive recent data favor a picture of a more disordered 
cluster which has not yet virialized. Slezak, Durret, & Gerbal (1994), Mohr, Fabricant, & 
Geller (1993), and McMillan, Kowalski, & Ulmer (1989) have found evidence of substructure 
from ROSAT and Einstein studies of the cluster. Andreon (1994) cites the non-uniform 
distributions of various Hubble types and along with Branduardi-Raymont et al. (1981), 
the observation that the optical and X-ray centers of A 426 do not coincide. The analysis 
of Schwarz et al. (1992) shows a cool subcluster in the eastern region of A 426, with an 
X-ray temperature > 2 keV less than the mean cluster average of 6-7 keV. This feature lies 
within 1 Mpc of NGC 1275. In other words, even though NGC 1275 is rather clearly sitting 
close to the dynamical center of its very rich surrounding region, that region has just as 
clearly not yet settled into equilibrium. If the Perseus cluster itself is an amalgamation of 
many subclusters, then we could reasonably suggest that NGC 1275 started out in one of 
these smaller clusters, in which lower- S'jv ellipticals are quite common (Harris 1991). But 
as a counterargument, examination of some prototypical high-Sjv galaxies such as M87 in 
Virgo (Binggeh & Tammann 1989), or NGC 4874 in Coma (Mohr et al. 1993) shows that 
they too belong to clusters which display plain evidence for substructure. 

These arguments lead only to possible links with galaxy environment and environmental 
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history. On the available evidence, we are forced to conclude that the specific frequency of 
a cD galaxy is not linked to the evolutionary state of its host cluster in any obvious way. In 
other words, Sn seems more likely to be a result of the initial conditions of the host gE 
galaxy. 

The most conservative synthesis of this discussion is to propose a scenario in which 
NGC 1275 originated as a large, or moderately large, E galaxy within a cluster that was 
significantly smaller than A 426 is today. There, proto-NGC 1275 would have had a GCS 
of average specific frequency. Over time, it may have merged with other galaxies, perhaps 
enough to cause a modest decrease in the number of GCs per unit luminosity. At the same 
time, it was fortunate enough to become the visible center of the larger Perseus cluster 
as other galaxy subclusters were amalgamated, building up to the system we now see. 
In this scenario, we would interpret the extensive nuclear activity in the center of NGC 
1275 as providing an intriguing window on cluster formation processes, but perhaps as 
inconsequential to the Sn of the system, since the globular clusters formed there will not 
survive long or add significantly to the old-halo population. 

This general interpretation probably makes the least demands on the ill-understood 
details of the merger process (either the merger rate itself or the cluster formation efficiency, 
for which we do not yet have a detailed physical theory). However, it is clear that other, 
more extreme, interpretations cannot be strongly ruled out. 

5. Summciry 

Our deep R-band survey of the old halo of NGC 1275 shows that this galaxy has an 
low-to-avcragc population of globular clusters despite its central location and dominant 
position within the Perseus cluster. Various arguments suggest that NGC 1275 might be 
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anything from a "starpile" of many previous mergers, to an initially enormous galaxy that 
has undergone only minor accretions since. In our view, the most likely possibility is that 
NGC 1275 started out with much the same specific frequency {Sn ~ 4) as it has today 
Results of this study indicate that GC formation is connected neither with a galaxy cluster's 
cooling flow, nor the evolutionary state of its cluster environment. Unfortunately, we are 
still left without any firm indications as to the formation mechanisms of high-S'Ar galaxies. 
Additional data for the NGC 1275 GC system in the form of accurate photometric color 
indices for the old-halo clusters would help to clarify some aspects of the argument. 
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Figure Captions 

Fig. 1. — The field surrounding NGC 1275 as imaged with CFHT and HRCam. The field of 
view is 2.7' on a side; this image is a composite of 20 300s exposures in the R band. 

Fig. 2. — Completeness of our i?-band photometry as a function of magnitude, as derived 
from artificial-star simulations and DAOPHOT. The two-parameter interpolation function 
(see text) is shown as the line through the data points. The limiting magnitude, defined as 
the magnitude at which / = 0.5, is at i? ~ 24.5. 

Fig. 3. — Image classification moments for the measured objects in the program field. Here 
Ri is an isophotal magnitude, ri a radial moment measuring image profile wingspread, and 
r_2 a radial moment measuring central peakedness (from Harris et al. 1991). Most detected 
objects in the frame are visibly nonstellar, with the 'starhke' objects showing up as the 
narrow sequence in the lower part of each panel. Objects with ri > 2.5 px are rejected from 
the sample; see text. 

Fig. 4. — Spatial locations of globular cluster candidates (faint, starlike objects) around 
NGC 1275. The elliptical annuh, centered on the galaxy, follow the shape of the galaxy 
light and are used to calculate the radial distribution shown in the next figure. Objects with 

22 < < 23 are represented by crosses, 23 < i? < 24 by filled triangles, and > 24 by 
squares. 

Fig. 5. — Radial distribution of globular cluster candidates around NGC 1275. Here a is 
the number of starlike images fainter than R — 22 per arcmin^, plotted as a function of 
radius (annular semimajor axis). The horizontal dotted lines show the adopted level of the 
background for a = —1.1 (dotted line), a = —1.3 (short dashes), and a = —1.5 (long 
dashes). The corresponding power law fits are very similar to each other. The innermost 
two annuli of the previous figure have been excluded from this graph due to poor number 
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statistics. 

Fig. 6. — The luminosity function (GCLF) for the old-halo globular cluster population in 
NGC 1275. The sohd symbols with error bars show the residual GCLF (j){R) (number 
per quarter-magnitude bin) after subtraction of background. The adopted background LF, 
defined by the three outermost annuli in Fig. 5, is shown for comparison by the dashed line. 

Fig. 7. — Radial distribution of the GCS, replotted for comparison with M87. The solid dots 
show the observed points as in Fig. 6 above. The labelled hue shows the distribution that 
would have been observed if NGC 1275 had the same specific frequency as M87, measured 
to the same limiting magnitude (data from McLaughlin et al. 1994). 
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TABLE 2: N1275 Globular Cluster Population^ 
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(1) Background subtracted and corrected for areal incompleteness. 



TABLE 3: Variation in Sj^ 



-1.1 (30.3 , 18.3) (25.1 , 21.8) 
-1.3 (33.4 , 14.9) (28.3 , 18.2) 
-1.5 (35.7 , 12.3) (30.7 , 15.5) 



380 ± 114 11 000 ± 3000 6.0 ± 1.7 
264 ± 80 7 800 ± 2100 4.3 ± 1.2 
190 ± 60 5 800 ± 1600 3.2 ± 0.9 



(a) unweighted values 

(b) weighted values 
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